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USER'S GUIP£ FOR SUBROUTINE PRNT3D 



Identification 

i.PRNT3D - A Subroutine Which Generates Two and Three Dimensional 
^ ' Printer Plots 

Author r Larry Gales 

Date July 1978, Center for Quantitative Science in Forestry, 
Fisheries and Wildlife, University of Washington 



^' Purpose- 

7 PRl'JTSD is a subroutii^e package which generates a variety of 

printer plot displays. The displays include single and multiple x vs. y 

I 

functions, multi-value x,y f^l^tionsbips, and density plots which 
simulate three-dimensional effects by means of overprinting. The 
package features one- and two-'dimensional interpolation, "zoom-in" 
capabilities, automatic scaling, logarithmic scaling, flexible titling, 
and multi-page plotting. Eacl^ plot fits on a standard 8-^ by 11 inch 
page with margins of sufficient size to permit incluj^iijii^in three ring 
bidders. Multi-page plots ^re automatically distributed over a number 
of -such 8% by 11-inch pages vitb sufficient annotation to permit recon- 
struction of the entire image. 

PRNT3D is built around an NX by NY grid of cells called an 

J^.^^lmage space. The image space is filled with numeric values from a 

binary file of x,y,z coordinates which contains data points as follows: 

X,, y,. z. 
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The X., coordinates of • each point specify its location in the x, y 
plane and the coordinate is interpreted as its height. Each x, y 
coordinate pair locates an appropriate cell within the NX by NY grid 
ipto which the z coordinate value is stored. Multiple z roapped 
to a single cell are either averaged together, or are re-, :r,- . .,d by 
the last encountered z value, according to an option se the 
user. 

Once the image space is complete, it is sent to an ouc^^Mt 
routine which printr it along with titles and appropriate scaling 
information. The association between a z value in the image c;race 
and its printed representation is determined as follows. Let 
ZICN < ZMAX be two values selected either by the user or by PRNT3D, 
and let Az = (ZMAX - ZMIN)/8. The z values are divided into 10 
levels as follows: 

Level 0: ^ 1 ZMIN 

Level 1: ZMIN < z < ZMIN + Az 

Level 2: ZMIN + Az < z < ZMIN + 2Az 

• • ' 

• • • 

Level 8: Z1«N + 7Az < z < ZliAX 

Level 9: ZMAX < z 

Each level is represented by a character or combination of 
overprinted characters. The user may opt to have the levels repre- 
sented directly by the single characters blank, 1, 2, 3, A, 5, 6, 7, 
8, or 9 orW the relative darkness of overprinted characters. The 
former is usually preferred for line plots while the latter is pre- 
ferred for surface plots. : 
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The following example should clarify the preceding discussion. 
Consider an image space constructed using the following parameters: 







D 








NY 




4 








XMIN 


B 


-100 


(minimum allowed 


X- 


-value) 


XMAX 


8 


400 


(maximum allowed 


X- 


-value) 


YMIN 




0 


(minimum allowed 


y 


-value) 


YMAX 




40 


(maximum allowed 


y- 


-value) 


ZMIN 




0 








ZMAX 




80 









This image contains 20 grid cells which are numbered from (1,1) at the 

lower left comer, which handles x and y values in the ranges: -100 < x <^0, 



0 < y < 10, to (5, A) at the upper 


right corner, which handles x and y values 


in the ranges: 300 < x < AOO, 30 


< 


y 


< AO. The Az value is (80-0) /8-10, 


and the z levels are: 












Level 


0: 






z 1 


0 


Level 


1: 


0 


< 


z < 


10 


Level 


2: 


10 


< 


z < 


20 


Level 


3: 


20 


< 


z <. 


30 


Level 


A: 


30 


< 


z < 


AO 


Level 


5: 


AO 


< 


z < 


50 


Level 


6: 


50 


< 


z < 


60 


Level 


7: 


60 


< 


z < 


70 


Level 


8: 


70 


< 


z <_ 


80 


Level 


9: 


80 


< 


z 




Consider the following file of x, 


y, : 


z 


coordinates : 


Sequence No. x-coord 


y 


•coord 


z-coord Grid Cell Mapped to 


1 -50 




25 






AS (1,3) 


2 150 




12 






75 (3,2) 


3 350 




32 






92 (5; A) 


'A -90 




25 






0 (i;3) 


5 500 




20 






36 



6 
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The image space, when filled from the above file, conceptxially appears as 
follows (the actual printed image lacks the full grid lines): 



AO 
30 
20 
10 
0 











9 


3 














8 

















-100 



100 



200 



300 



AOO 



Note that the data file contains five points, but the inaga space shows 
only three. This is because (a) PRNT3D ignored point 5 since its x coordinate 
exceeded XMAX, and (b) by default, it averaged the z values (A5 and 0) for 
points 1 and A since they were mapped to the same grid cell, (1,3). The 
average value, 2=22.5, is represented by z level 3, and the other z values, 
75 and 92, are represented by z levels 8 and 9, respectively. If the image 
were displayed using overprinting, it would appear as follows: 



AO 
30 
20 
10 
0 



-100 



100 



200 



300 



AOO 



where grid cell (3,2) contains "W", "T", and "-" overprinted, and" (^5, A) 
contains "B", "M", "*". 
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PRNT3D contains several options which partially overcome the 
limitations of low resolution, sparse data representation, and device 
dependent printing characteristics, which commonly plague printer plot 
displays. 

The problem of low resolution is handled by the zoom-in and 
multi-page options. Zoom-in permits a user to blow-up a region of the 
data into as much detail as is desired by specifying the x, y coordi- 
nates of a rectangular window which encloses it. The data in this 
region are mapped onto the full NX by NY grid - data outside the region 
are ignored. Normally, a user first views the data set as a whole and 
then selects windows where more detail is desired. 

The multi-page option permits an effective resolution of up 
to a million picture elements over all or part of the data. Multi- 
paging is controlled by the parameters NX and NY which specify the 
number of x and y axis cells in the image space. If NX < 60 and 
NY < 45 the plot is restricted to a single page, otherwise it is 
automatically distributed over a number of pages. The maximum reso- 
lution is NX = 999, NY = 999. 

Sparse data presents a problem in printer plot displays, 
because it is difficult for a user to perceive the structure of ob- 
jects which are represented by a few scattered points. PRNT3D con- 
tains one- and two-dimensional interpolation or point enrichment 
options which add equi-spaced, linearly interpolated data points to 

the points in the original data file. Both types of enrichment -require 
that data sets be organized in special ways. The requirements an^ de- 
fails of enrichment are fully described in the appendix. 

er|c ^ 
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The problem of device dependence affects the visual quality of 
density (overprinted) plots because print combinations which are dark 
on one device may be relatively light on another. PRNT3D permits 
the user to map any of the 10 z value levels occurring in the data on 
to any of 10 predefined print combinations by means of an array named 
ZMAP. If ZMAP(i) = j, where 1 < i < 10 and 0< j < 9, then the ith z 
level will be represented by the print combination normally reserved 
for level j. For example, if the third print combination on one device 
were actually darker than the seventh, the user could interchange them 
by setting ZMAP(3) = 7 and ZMAP(7) = 3. ZMAP can also be used to reduce 
the number of printed levels and hence sharpen the contrast, to reverse 
the representation of levels (e.g. by setting Z\S^ « 9, 8, 7, 6, 5, 4, 
3, 2, 1, 0), or to display rough contour levels (e.g. by setting ZMAP = 
0, 1, 0, 1, 0, 1, 0, 1, 0, 1). 
Input 

This write-up assumes that PRNT3D is already incorporated as 
a subroutine in an existing program which uses the free-form input package 
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(Gales and Anderson, 1978) and follows the conventions outlined in the 
Design Standards for Computer Programs (Gales, 1977). It discusses 
only those plot directives which are assigned by the user through the 
free-form package. For a more thorough discussion of PR1^T3D refer to 
its programmer's guide. 

The user assignable variables in PRNT3D are XMIN, XMAX, YMIN, 
YMAX, ZMIN, ZMAX, XRICH, YRICH, DFAULT, OVPRNT, AVE, INT2D, NX, NY, 
and ZMAP, and are described as follows: 

^ INPUT TABLE 



NAME 



TYPE AND DIMENSIONS 



RANGE LIMITS 



DESCRIPTION 



XMIN 
XMAX 
YMIN 
YMAX 



ZMIN 
ZMAX 



Real 



Real 



^29 



29^ 



(-10"', 10"') 



29 29 
(-10^^, 10^^) 



XMIN, XMAX and YMIN, 
YMAX define a rectangu- 
lar window which encloses 
the data in the image space 
to be displayed. Data 
outside the window are 
not shown. If XMIN>XMAX 
and/or YMIN>YMAX the 
computer ignores them 
and constructs a window 
which just encloses all 
data in the data file. 

ZMIN and ZMAX are the 
lower and upper bounds 
for 10 levels which 
determine the printed 
representation of z 
coordinate values in 
the binary data file. 
If ZmnzZmX, the com- 
puter Ignores them and 
assigns the lowest and 
highes z values in the 
data file to ZMIN and 
ZMAX, respectively. 
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NAME 



TYPE AND DIMENSIONS 



XRICH 
YRICH 



Real 



DFAULT 



Real 



OVPRI^T 



Logical 



AVE 



Logical 



RANGE LIMITS 



(0, 10^^) 



(-10^^, 10^^) 



.T. or .F. 



.T. or .F. 



DESCRIPTION 



XRICH and YRICH control the 
Ax and Ay Increments 
used In both one- and 
two-dimensional inter- 
polation (enrichment) • 
If XRICH « 0 and/or 
YRICH « 0, no enrich- 
ment takes place. The 
user should note that 
the values of XRICH 
and YRICH should be 
coordinated not only 
with the data, but with 
the size of the window 
set by XMIN, XMAX, YMIN, 
YMAX. If XRICH or YRICH 
are too small, the en- 
richment process will 
consume too much compu- 
ter time, whereas large 
values of XRICH, YRICH 
will leave gaps. 

DFAULT is the default 
value assigned to all 
cells in the image 
space. DFAULT is 
usually set to zero. 

If OVPRNT is true, the 
z values in each cell 
in the Image space will be 
represented by a set of 
overprinted characters, so 
that high z values will 
appear dark(the lowest level is 
always blank). If OVPRNT 
is false, then the z values 
will be represented by one 
of the characters blank, 1, 
2, 3, 4, 5, 6, 7, 8, or 9. 

If AVE is true, then all z 
values mapped to a single 
cell in the ima^e space will 
be averaged. If AVE is 
false, the last z value 
stored in the deil takes 
effect. 
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NAME 



INT2D 



NX, 
NY 



ZMAP 



ERIC 



TYPE DIMENSIONS 



Logical 



RAl^GE LIMITS 



.T. or .F. 



Integer 



Integer (10) 



2 ^ 
2 < 



|NXl 
iNYl 



999 
999 



DESCRIPTION 



0, 9 



If INT2D is true, then 
two-dimensional inter- 
polation or enrichment 
will be applied to the 
binary data file, pro- 
vided that both XRICH 
and YRICH are greater 
than zero. If INT 2D 
is false, then one- 
dimensional interpola- 
tion will be applied 
if both XRICH and 
YRICH are greater than 
zero. Note that inter- 
polation can only be 
applied to data which 
are correctly organ- 
ized on the binary 
data file. 

NX and NY are the number 
of X and y cells in the 
image space. If NX<60 
and NY<45, the image is 
printed on one page, other- 
wise it is automatically 
spread over a number of 
pages. For multi-page' plots 
NX should be an exact multi- 
ple of 60 and NY an exact 
multiple of 45. If NX (NY) 

0, the number of cells is 
computed as the absolute value 
of NX(NY) and the •x(y)-axis 
is scaled logarithmically as 
are the points in the data set« 



The ZMAP array maps a 
given z level into one 
of the 10 predefined 
print combinations. 
Normally, the array 
ZMAP«0,1, 2,3,4,5, 6,7,8,9 
which means that the 
lowest z level is repre- 
sented by a blank, the 
next z level either by 
a "1" or a the 
next z level by a '-2" 



the 



or an - , 
last z level by a 9 or 
the overprinted set "B'*, 
"M", and 
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NAME 


TYPE AND DIMENSIONS 


RANGE LIMITS 


DESCRIPTION 


ZMAP 
(cont) 


Integer (10) 


0, 9 


However, if the user 
sets ZMAP(i) = j for 
any i or j then the 
1th z level will be 
represented by the jth 
print combination. 



Output 

The output for PRNT3D consists of the printed representation 

of the image space, four lines of text at the top of the plot which 

serve as a title, x and y axis text labels, x and y numeric labels of 

the form + N.NNN, and a legend underneath the plot which specifies the 

scale factors for the x, y, and z axes, whether x and y scales are 

linear (default) or logarithmic, the values for each of the 10 z-axis 

levels, and the number of points mapped to each z level. Since the 

latter figures are restricted to two digits, values in excess of 99 

are represented by -9. For example, if the plot legend were to read 

SCALE FACTORS » X-AXIS :E+06, Y-AXIS (LOG) :E-03, Z-AXIS:E+04 
Z0-Z4 = 3.219(17), 3.319(11), ... 
Z5-Z9 = 3.719(92), ... 

then all x-axis numeric labels must be multiplied by 10^, all y-axis 

values by 10"\ and all z-level values by lO"*; 17 points (whose z-axis 

values were < 32190) were mapped to level 0 (blanks), 11 points (whose 

z axis values were in the interval (32190, 33190]) were mapped to z level 

2, and so on. 

As an example, consider the first sample run whose image 
space contains 8 digits, 1 through 8, scattered sparsely over the plot. 
In this case, the four title lines and the x and y axis text labels are 
absent. Consider the point "8" which is located at the bottom nd^dle 
of the plot. Its X and y axis coordinates are (x > 2.996, y = 1.000) 
and since the x and y axis scale factors are E-fOO (corresponding to 
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10® » 1) the true coordinates remain (x«3, y=l). More than 99 points 
were mapped to z level 0 (which always show up as blanks), 1 point each 
was mapped to z levels 1 through 8, and none were mapped to z level 9. 
All scales were linear, so the word "LOG" did not appear in its legend. 
Restrictions 

The one- and two-dimensional enrichment functions should only 
be used on data which are correctly organized on the binary file of 
X, y, z coordinates. 
Error Messages 

PRNT3D issues four types of fatal error messages: range check 
errors, errors arising from enrichment of faulty input data, errors due 
to zero width windows, and errors arising from non-positive data values 
when logarithmic scaling is used. The first type are of the form: 

ERROR NO. X IN PRNT3D 

yy. . .y OUT OF RANGE 
yy. . .y « dd. . .d 

where x is 1, 2, or 3, yy...y is the name of a variable (NX, NY, or ZMAP) 
and dd...d is the value of the variable. The second type are of the form: 

ERROR NO. 4 IN PRNT3D 

A TRIANGULAR REGION IS MALFORMED 
ITS POINTS ARE: , 
(Xj^, y^^, z^), (x2, y2» ^2^* ^^3* ^3* ^3^ 

ERROR NO. 5 in PRNT3D 

END OF FILE ENCOUNTERED BEFORE • \ 

A TRIANGLE IS COMPLETE 



Error number 4 occurs if the three vertex points of a triangle 
form a single line, and number 5 occurs if the last triangle possessed 
only one or two vertex points. The third type of messages are of the form: 

ERROR NO. X IN PRNT3D 

COMPUTED a-AXIS WINDOW IS ZERO WIDTH 
AT a = dd . . . d 

where x is 6, 7, or 8, and a is X, Y, or Z. These errors occur if the 
user is letting PRNT3D determine the size of the window (by setting 
XMIN > XMAX and/or YMIN > YMAX) and all the x coordinates are equal or 
all the y coordinates are equal (for example, if only one point appears 
in the image space). The last type of error messages are of the form: 

ERROR NO. X IN PRNT3D 

LOG OF NON POSITIVE a VALUE 

AT cxMIN « dd...d 
where x is 9 or 10, and a is X or Y. 
Sample Runs 

The annotated listing on the next few pages illustrates the 
control cards and input cards for several sample runs. The sample runs 
show the effects of various plot options applied to a single set of x, 
y, z coordinates on a binary data file. 
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•ACeDUNT#9XMO8C00# . 

' COMMENT. ************************************* 
COMMENT. ♦ THE FIRST CARD ABOVE IDENTIFIES * 
COMMENT. ♦ THE JOB* SPECIFIES THE MEMORY * 
COMMENT. ♦ REQUIREMENTS (55000 OCTAL) AND THE* 
COMMENT. ♦ ESTIMATED CENTRAL PROCESSOR TIME * 
COMMENT. ♦ IN SECONDS (50 OCTAL). THE SECOND* 
COMMENT. * CARD IDENTIFIES THE BUDGET AND * 
COMMENT. <* PASSWORD. * 
COMMENT. ************************************* 

COMMENT. 

A,TTAC'H#BPR30#I0«BPR30. 

ATTACH#BFF#ID«BFF. 

MNF#L*0#E«1#B«LARGE. 

COMMENT . ************************************* 

COMMENT. ♦ BPR3D IS THE FILE CONTAINING THE ♦ 
COMMENT. * 3D PRINTER PLOT SUBROUTINE IN * 
COMMENT. ♦ BINARY FORM» BFF IS THE FORMAT * 
COMMENT. ♦ FREE INPUT SUBROUTINE IN BINARY * 
COMMENT. * FORM. THE MNF CARD COMPILES THE * 
COMMENT. ♦ PRNT3D TEST PROGRAM AND WRITES * 
COMMENT. ♦ THE BINARY PROGRAM TO FILE LARGE. * 
COMMENT. ************************************* 
COMMENT. 

LOA0#LARGE« BPR3D#BFF. 
EXECUTE. 

COMMENT. ************************************* 
COMMENT. ♦ THE LOAD CARD LOADS THE TEST * 
COMMENT. * PROGRAM (ON FILE LARGE) INTO MEM- * 
COMMENT. * ORY# ALONG WITH BPR3D AND BFF. ♦ 
COMMENT. ♦ THEN CONTROL IS PASSED TO THE TEST* 
COMMENT. * PROGRAM* WHICH BEGINS EXECUTION. * 
COMMENT. ♦ INPUT IS HANDLED BY BFF# AND * 
COMMENT. * OUTPUT IS PLOTTED BY BPR3D. * 
COMMENT. ************************************* 
♦EOR 



♦♦***PROGRAM GOES HERE***** 



*EOR 

/***********t^***************** PUN I **************^^^t**************** 



I THE FIRST RUN PLOTS EIGHT POINTS IN THE REGION I < X < 5 AND 

/ 1 < Y < ^. THE POINTS ARE NUMBERED 1 THROUGH 8 AND NO 

/ INTERPOLATION IS DONE BETWEEN POINTS. THE DEFAULT X# Y# AND Z 

/ COORDINATES FOR THE 8 POINTS ARE AS FOLLOWS 

/ POINT X 

/II 

/ 2 3 

/ 3 1 

/ A 3 

/ 5 4.2 

O 
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3 
2 

2.4 
2.8 



1 
1 
1 

1 
1 



/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 

/ 
/ 
/ 



6 
7 

6 



5 
3 



3 1 
2.5 1 
1 1 



THE POINTS ARE DESCRIBED BY A TWO DIMENSIONAL ARRAY 

NAMED -COORD- WHICH IS STRUCTURED AS FOLLOWSl 

C00RD(1#P) ■ X COORDINATE FOR POINT P 

CaORD(2#P) ■ Y COORDINATE FOR POINT P 

Ct>0RD(3*P) ■ Z COORDINATE FOR POINT P 

FOR EXAMPLE* SETTING C00RD(2*1) " 

Y COORDINATE FOR POINT 1 EQUAL TO 

C0bRDtl*5) • <».2 SETS THE X COORDINATE OF POINT 

5 EQUAL TO <».2. 



3 SETS THE 
3, WHEREAS SETTING 



IN ORDER TO NUMBER THE POINTS IN THE PLOT* 

THE DEFAULT Z COORDINATES FOR ALL POINTS* THAT IS* 

C00RD(3»J) FOR J ■ 1* ...» 6* APE CHANGED FROM 1 

TO THE VALUES 1 THROUGH 8* AS FOLLOWS 



C00PD(3*1) - 1* 
C00RD(3*^) ■ <»* 
C00RD(3*7) ■ 7* 



CQ0RD(3*2) ■ 2* 
C0DRD(3*5) ■ 5* 
C00RD(3*8) ■ B» 



C00RD(3*3) - 3* 
C00RD(3*6) - 6* 



THE FOLLOWING PLOT PARAMETERS ARE INPUT 

YRI CH 



XRICH ' 
INT2D ' 
OVPRNT 



0* 
.F.* 
► . F .* 



0* 



THE SECOND RUN DOES A ONE-DIMENSION INTERPOLATION BETWEEN 
SELECTED POINTS FORMING 3 TRIANGLES IN THE REGION. THE TRIANGLES 
FORMED HAVE THE FOLLOWING POINTS AS VERTICES 



TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 



II 
21 
31 
41 
51 
6 I 
71 
81 



POINTS 
POINTS 
POINTS 
POINTS 
POINTS 
POINTS 
POINTS 
POINTS 



1*3* 
2*3* 
2*<»* 
2*5* 
5*6* 
5*8* 
5*<»* 
<**3* 



AND 
AND 
AND 
AND 
AND 
AND 
AND 8 
AND 8 



/ THE DEFAULT X* Y* AND Z COORDINATES FDR THE POINTS ARE THE SAME 
/ AS FOR RUN Ic THE FOLLOWING PLOT PARAMETERS ARE INPUT 

XRICH ■ 0.05* YRICH ■ 0.05* 

INT2D ■ .F.* 

OVPRNT ■ .F.* $ 
/ ' 
/ NOTE! THE INTERPOLATION ALGORITHM DEPENDS ON THE ORDERING OF THE 
/ POINTS IN THE DATA SET. THE ORDERING OF THE POINTS IN TrtE DATA SET 
/ FDR THIS PARTICULAR RUN IS SUCH THAT NOT ALL TRI ANGLES . AllE 
/ COMPLETED BY THE INTERPOLATION. THUS THE LINES BETWEEN POINTS I 
/ AND 2 AND BETWEEN 2 AND 6 DO NOT APPEAR IN THE PLOT. 



17 



r% -15* 

/tTHE THIR-£^^ RUN IS A DETAIL Of RUn 2 AND DISPLAYS THE REGION 
/ DEFINED BY 3 < X < 5 AND I < Y * 3 
/ THE FOLLOWING PLOT PARAMETERS ARE INPUT 
/ - 

INT2D .F.# 
^ OVPRNT ■ .F., 
^ XHIN • 3* XMAX ■ 5f 
•V^N ■ 1, YMAX -3, 
-rj^!9H„S« 0.01* YRICH - 0*01, I 

I 

I THE FOURTH RUN DISPLAYS A SURFACE WITH INTERPOLATION. THE 
/ Z COORDXNATES FOR POINTS Z» k, AND 5 ARE RAISED TO A VALUE 
/ OF 9. THE .REMAINING POINTS HAVE THE DEFAULT Z COORDINATE 
/ VALUE OF' T. 
/ 

COORDC.3,2) ■ 9, C00RD(3^<r) - 9, C00RD(3»5) ■ 9, 

/ 

/ THE FOLLOWING PLOT PARAMETERS ARE INPUT 
/ 

XRICH ■ 0.035, YRICH * 0.035* 1NT2D " .T., 
OVPRNT » .T., $ 

/ 

I 

I THE FIFTH RUN IS THE SAME -AS THE FOURTH RUN WITH REVERSED VIDEO. 
/ 

C00RD(3,2) ■ 9, C0QRD(3/<r) • 9, C00RD(3,5J ■ 9, 

/ 

/ THE FOLLOWING PLOT PARAMETERS ArE INPUT — 
/ 

OVPRNT ■ .T., 

XRICH 0.035, YRICH * 0.035* InT2D ■ .T,, 
ZMAP ■ 9,8,7,6,5,^,3,2>1»0, » 

/ 

/*«4i*4,4i:««««««««««*«<»« «««*4>4. *#% STOP PROGRAM «««««««4>««*«««4>«***«*««««>#* 
/ 

F-INIS ■ .T., $ 

/ 

♦ EOR • 

♦ EOF 



\ r 



18 



PROGRAM -TEST- READY FOR INPUT 



f ^^^^^^^^^^^^^t*************** RUN 1 ********************************* 



I 
I 
I 
I 
I 
I 
I 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 



/ 

/ 



THE FIRST RUN PLOTS EIGHT POINTS IN 
1 < Y < THE POINTS ARE NUMBERED 

INTESPOLATION IS DONE 
COORDINATES FOR THE 8 



THE REGION I < X < 5 ANH 

1 THROUGH 8 AND NO 

BETWEEN POINTS. THE DEFAULT X» Y» AND Z 
POINTS ARE AS FOLLOWS 



POINT 

1 
Z 
3 
^ 

5 
6 
7 
8 



1 
3 
1 
3 

5 

3 



3 
<» 
2 

2, 
2 
3 

2 .5 

1 



,9 



1 
1 
1 
1 
1 

1 



THE 



.... POINTS ARE 
NAMED -COORD- 
COORDdf P) 
C00RD(2»P ) 
C00RD(3»P) 
FOR EXAMPLE* 
Y COORDINATE 
C00RD(1»5) • 
5 EQUAL TO <» 



DESCRIBED BY A TWO 
WHICH IS STRUCTURED 
■ X COORDINATE FOR 
• Y COORDINATE FOR 
. ■ Z COORDINATE FOR 
SETTI^4G C00RD(2*1) ■ 
FOR POINT 1 EQUAL TO 



4.2 

2 



DIMENSIONAL ARRAY 
AS FOLLOWS! 
POINT P 
POINT P 
POINT P 
3 SETS THE 
3* WHEREAS SETTING 



SETS THE X COORDINATE OF POINT 



IN ORDER TO NUMBE« THE_POINTS 
THE DEFAULT 
C00RD(3>J) FOR 
TU THE VALUES 



COORDINATES FOR 
J ■ \f ...» 8> 

1 THROUGH 8* AS 



IN THE PLOT* 

ALL POINTS* THAT IS* 

ARE CHANGED FROM 1 

FOLLOWS 



C0aRD(3*l) • 1* 
C0nRD{3*4) ■ A* 
C00RD(3»7) ■ 7, 



CO0RD( 3* 
C00RD(3* 



C00RD(3* 8) 



2* 
5* 

8* 



C00RD(3* 3) 
C00RD(3*6) 



3* 
6* 



THE FOLLOWING PLOT PARAMETERS ARE INPUT 



XRICH ' 
INT2D ' 
OVPPNT 



0* 
.F 



YRICH ■ 0* 
$ 
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1.000 
X— 

A. 000' 



3.659 



1.610 

• — X — 



2.288 
, — X — 



3.318Y 

2.977 

2.636 

2.295 

1.955 

l.Sl'i 

1.2 73 

1.000 
1. 



SCALE F 
Z0-Z4- 
I Z5-Z9- 



2.966 
X--' 

2 



3.6AA 
X 



4.322 

X — 



5.000 
• — X 



X-- 

00 



X 

1.610 



2.2S8 



9 

2.966 



3.6A'. 



y— . 

A. 322 



• — X 
5 .000 



CTGRS • X-AXIS: E+00 

0(-9), 1.000 { L)f 
5.000{ 1>* 6.000{ l)r 



Y-AXISJ E+00 Z-AXISI 
2.000( 1)* 3.000{ 1)» 
7.000( D* 8.000( 1 )* 



E+00 

A. 000 { 
9.000( 



1) 
01 



ROGRAM -TEST- READY FDR INPUT 



RUN 2 ♦♦♦♦♦♦♦♦♦**♦♦*♦♦**♦♦♦***♦*♦*♦♦*** 



THE SECOND RUN DOES A DNE-DI MENS I DN I NTEP POL AT 1 DN BETWEEN 
SELECTED POINTS FORMING 8 TRIANGLES IN THE REGION. THE TRIANGLES 
FORMED HAVE THE FOLLOWING POINTS AS VERTICES 



TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 
TRIANGLE 



11 
21 
3t 
<»: 
5t 
6t 
7i 
81 



POINTS 
POINTS 
POINTS 
POINTS 
POINTS 
POINTS 
PO INTS 
POINTS 



1*3, 
2,3, 
2, A, 
2,5, 
5,6, 
5,8, 
5, A, 
A, 3, 



AND 
AND 
AND 
AND 
AND 
AND 
AND 8 
AND 8 



ERIC 



THE DEFAULT X, Y, AND Z COORDINATES FOR THF. POINTS ARE THE SAME 

20 



-AO- 

J AS FOR RUN I. THE FOLLOWING PLOT PAR Aht , . S ARE lN«*Uf 

XRICH ■ 0.05» YRICH ■ 0.05* 

INT2D • .F.* 

OVPPNT ■ .F., S 



■ 1' 

O 

ERIC 



-19- 



1«000 
X— 

4.000' 



3.659 



1.610 



2.238 



2.966 



3.318 



2.977 



2.636 



2.295 



1.955 



1.61A 



1.273 



1.000 



1 
11 
1 



11 
1 



1 

1 
11 



1 
1 



1 

111 
111 
1 1 
I 
I 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3.6AA 
X-- 



<».322 
. — X— 



5.000 
• — X 



11 



11 



1 
1 



11 



1 

11 



11 
1 



11 



1 

1 



1 
1 



1 

11 



il 
1 



11 

1 



1 
1 



11 



1 

1 
11 

11 111 
11 1 111 
1 1 111 
111 
1 

11 

1 

. ^1 

Si 



1111 
111 1 
1 
1 
1 
1 
1 
1 
1 



1 
1 
1 
1 
1 

1 1111 
111 



111 



111' 



1 1111 
1111 
111 1 11 



1 
11 
1 



l' 

1' 



11 
1 



1 
1 

111 
111 
IIH 

1 



1 

11 
1 

1 

11 

1 11 
1 1 
1 11 
1 1 
1 1 
1 11 
11 
111 



11 
1 



11 



11 
1 



1111 
11111 1 
11 
1 

11 
1 1 
1111 
1 

11 



X— 
1.000 



--X — 
1.610 



2.28'j 



2.966 



X— 

3.6*1^ 



X— 

<».322 



•X 

00 C 



SCALE FACTORS ■ X-A>.IS» E+00 
ZO-Z«»- 0(-9), 1.000 (-9)* 

Z5-Z9* 5.000( 0 )* 6.0G0( 0)* 

ROGRAM -TEST- READY FOR INPUT 



Y-AXISI E+00 Z-AXISl E+00 
2.000( 0)* 3.000( 0)» <».000( 0) 
7,000( 0)* 8.000( 0), 9.000( 0) 



NDTEt THE INTERPOLATION ALGORITHM DEPENDS ON THE ORDERING OF THE_ 
POINTS IN THE DATA SET. THE ORDERING OF THE POINTS IN THE DATA SET 
FOR THIS PARTICULAR RUN IS SUCH THAT NOT ALL TRIANGLES ARE 
COHPLETED BY THP INTERPOLATION. THUS THE LINES BETWEEN POINTS 1 
AND 2 AND BETWEEN 2 AND 6 DO NOT APPEAR IN THE PLOT. 

RUN 3 ♦♦♦***♦*****♦»♦**♦♦♦:♦*♦*♦♦♦»♦****♦ 



/ THE THIRD RUN IS A DETAIL OF RUN 2 AND DISPLAYS THE REGION 
DEFINED BY 3 < X < 5 AND 1 < Y < 3 ; . 

THE FOLLOWING PLOT PARAMETERS ARE INPUT 



/ 
/ 
/ 



INT2D 
OVPRNT 
XMIN ■ 



» . F • » 
■ • F . » 

3, XMAX 



;6i?H-.^&.oi:"*S5i?fi - 0.01, s 



3.000 
X- 
3.000Y1 
II 
!1 
II 
II 

2.773Y1 

ii 

II 
II 

2.545Y1 

II 
II 

1 1111 
II 

2.318Y1 
II 

11 

U 

2.091Y1 
II 
II 

\l 

l.e6<rYl 
II 

u 
II 
II 

1.636Y1 
II 
II 
II 

II 

1.^09Y1 

H 

II 
II 

1.182Y1 
1 1 

mil 
nil 

l.OOOYll 

X 

3.000 



3.305 
X 



3.64<> 
• — X 



3.983 

— X 

11 
11 



<i.322 
X — 



A. 661 
X 



11 
11 



mill 



nil 
mil 
11111 
mil 
11111 
mil 



1 

11 
11 
11 



111 

11111 ^1 

11 
11 
11 
11 
11 
1 



11 
11 
11 
11 
11 

• 11 1 
11 11 
1 11 
1 111 
1 11 
11 11 
11 1 
nil 
nil 
nil 



1 

1 1 



11 
11 
11 



11 
11 
1 



11 
11 
11 



11 
11 
n 



11 
11 
11 



11 
11 
11 



11 
11 
n 



5.000 
X 

nil 
11111111 
mm 11 
111 
11 
11 

il^ - 
11 11 
11 n 
11 11 
11 

11 
11 
11 



— X — 

3.305 



3.6<i<i 



■ X— 

983 



X— 

A. 322 



— X 

4.661 



■X 

,000 



JCALE FACTORS ." X-AXISLEtOO YrAXISi 



25-Z9- 



0(-9)* 
5.000? 0)* 



1.000 (-9), 
6.000( 0)t 



ROGRAM -TEST- READY FOR INPUT 



_^ E+00 
2.005( 0)* 
7.000( 0)* 



Z-AXISi E+00 
3.000( 0 )» 



8.000( 0)* 



A.000( 0) 
9.000( 0) 



f 
I 
I 
I 
I 
I 

I 
I 
I 



RUN A 

THE FOURTH RUN DISPLAYS A SURFACE WITH INTERPOLATION. THE 
Z COORDINATES FOR POINTS 2, A, AND 5 ARE RAISED^TO A VALUE 
OF 9. THE REMAINING POINTS HAVE THE DEFAULT Z COORDINATE 
VALUE OF 1. 



C00RD(3*2) ■ 9# C00RD(3*<») ■ 9# COORD(3*5) 
THE FOLLOWING PLOT PARAMETERS ARE INPUT 

0.035# INT2D ■ .T., 



9, 



XRICH « 
GVPRNT 



0.035, 
• . T. > 



YRICH 
S 



ERIC 
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1.000 
X — 

A.OOOY 

r 



3.659 



3.318Y 



2.977Y 



2.636Y 



1.95? 



1.61^Y 



1.273 



l.OOOY 

X — 
1. UOO 



«- ■ 0 ■ ■ • 



1.610 2.288 2.966 3.64't A. 322 5.000 

^ ■ 

iff f i 

UWUVif BBIiiWU 
^UUWUUtilBBBf tWU# 

AA^M###UUUIfBBIBBBtVUU##A 
A&AAA#^###UWiiB8BBBBBifUU#AAA 
**k AAAAA####UUif BBBB8IBBffUU#AAA« 
AAAAAA###UUffBBBBHBBBB{f UW^AAA** 
. AAA AAA##UUIfBB6BBeBBBBi«UUM A A«^. 

«***««AAAAA##UUIBBBBBBBBBBBtifU#^#AA«. . . 

^♦^♦♦♦AAAA^^UUliBBBBBBBBBBBifffWU^^AA*. .-- 
I I! .«««««*AAAA#VU{iBBBBBBBBBBBBiVeuu##AA*. 

♦♦♦♦♦♦AAA^UU«iBBBBBBBBBBBBBii»UU##AA*. 

A A^UUViBBBBRBBBBBBBBBif iUU^#AA«. . - > — 

.a ««4-«AA#WUf iBB6BRBBBBBBBBBBBilUV##A**.- — 

..- ♦♦*AA#UWilBBBBBBBBBiBBBBBBBllU##A*.«i- 

AA#UUf iBBBBBBBBBBBBRBBBf Bi##** .-- 

««AA^UUff»BBBBBBBBBBE8B<9ifiuyAA«.-~ 

.««AA#UUifBflBBBBBBBBtfffiieuU#A«.-- 
-....• ■ . . . ««AA#yUiiBBBBBflBffiiitiSiiU#A*. - — 

.«*AA^UUifB8BBltiiiiffiUUU^^&*>-- 

A A #WyiIBtlf Iff ««UUUWU##A«.- — 

««AA#UUf fifiVVVUUUUUUVUAA*. 

■-. .♦♦AA#UWUHiiWWUUWUUU##AA*."- 

■-. , ♦♦AA^^WWUUWWUWWWUM^^AA* . ■ — 

... .♦♦AAA#^UWWUWWUM#Mf A* .■- 

... . .♦♦AAA#3UWU*#?##^V^AA*. ■- 

... ..♦♦AAA^^^^^^^^^AAA*. .■- 

— .....♦♦♦AA######AAAAA**."- 

... ..♦♦AAAllAAAAAAAA*. ■ — 

— .....♦♦♦AAAAAAA 

AAA****** .■- 

■ ■■. .♦♦*♦♦♦****«■ — 

.*♦♦***. 
. ■ — 



— X 

1.610 



2.288 



2.966 



3.6AA 



. — y 

'».322 



• — X 
5.000 



!S^^z^/^^^°^^(-9?;^M5io§r-§), 3.Uf^!,^Tooo(-9) 

25-29- 5.0001-9)' 6.000{-9), 7.000(-9), 8.000{-9), 9.000(-9) 

PROGRAM -TEST- READY FOR INPUT 



RUN 5 



^♦♦♦♦♦♦♦♦♦♦♦♦♦***** ♦+**♦****♦* 

/ THE FIFTH RUN IS THE SAMF AS THE FOURTH RUN WITH REVERSED VIDEO. 

' C00RD(3»2) • 9» C00RD(3»<.) • 9, C00RD(3*5) ■ 9, 

/ THE FOLLOWING PLOT PARAMETERS ARE INPUT — 

/ ^ \ 

XPICH^-'o:035, YRICH - 0.035, ^ INT2D - .T., 
ZMAP ■ 9, 6, 7, 6, 5, 4, 3,2,1,0, S 



or 



ERIC 



I'Ooo 1:$12 --X-- ^*X^^ ^'P^ ^'y^° 



iBB8eBeiBB8Bf!eefleflB8BeRBf!BBeil 
■iBBBBBBIBBBBBBBBBRBBBBBBRRK 



— — Da,, 



1<A 



.'se.ii?*;':^ 



^.OOOYBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
IBBBBBBBBBBBBBBRBBBIBBBBBBBBB — 

IBB8BBBBBBBBBBBBBSBBBBBBBBB 

liBBISniiBiilBBBililf RBBa•>■-<— 
IBflBBBBBBBBBBBB6BBBB8BB.----•— ' 

3.659YBB8BBBBBBBHBBBBBBBBB — 

IBBBBBBBBBIBBBBBBBB^* — 

IBflBBB8BBBBBBBBBB+«««« — 

IBBBBBBBBBBBBBB A A« . . .-- — 
I BBBBBBBBBBBBAAAAA* .•- — 

3.318YBBBBBBB8BB^#AAAAA A««*«««. .-- — 
IBBBBBBBB#^^##AAAAAA«««««. .-- — 
IBBBBBBUUU###y#AAAAAA««*«. .-- — 
IBBflBiUUUUU#^###AA AAAA«*«« .-- — 
IBBiffiWUUUUU#####AAAAAA«««. - — - 

2.9 77YlfffffiiWUUUUW^^##^#AAAAA««.-- — 
lilSllff UUUUUW###^^^AAAA«« .--— 
IBBBiifiCUWUWUW##^##^AAA««. — 
IBBBBtf iiiWUUUUW#4M##AAA««.- — 
IB8BBBif ViffWUWUW#^ ###^AA«« .-- — 

2.63&Y8BBBBBf iiiiWUUWUW#^#i'AA«« .>- — 
IBBBBBBBffifffUWWWWVl###AA«« .-- — 
IBBBBBBBBif iiiWUUWUU##AA««. — 

IBBBBBBBBBiliifJWWUW^^AA**."" ■■■■"li'*ii'.SS5J:'S2siii. 

IBBBBBBBBBBiifSitfW WW#^AA** .■■ •**A#WiB091BBB~ 

2.29 5YBB88BBBBRBBfi»»ifJWW##AA**."-- • 1* ♦^.^1)(^2I2S99§H;? 

IBiBBBBBB8BBBV«V«SUU##AA«« . A^WSGSgB8:33B'^F. ^r- 

lB8BBBflBBflBBB8ffiiyw##AA***. 5*^''i(f 55i52iSiviP-.r [-^ 

IBflBBBBBBBBBBRRiIiWU#^^AA*«^. . « « A ^WV8dB8?fleB9S?i^'- 

IBBBBBBBBBB88BBBfffffUU###AA«*« «««A##W«BBBBBBE;8Bt)SR^' ? 

1.953YBeiiBBBBBBBBSRBBiiWUU##AAA«« «««««AA#uieBBBBBBB8B36RS-.^ 

iSSSUUSss'.ssssyssiii^.ryiSirArAiir/AiSiSHSissussisuiii-d^ 

1 f,l<.YB8B8BBB8BBBBBBBBIBBBRBVUUU##AAAAAA##UkiBBBBBBiBBBBBBIBB6nB^3 
IBBBBBB8BBBBBBBBBByBflBffWWUi>#|iA###|l|>WfSBflBBBBBBBBBBBB^ 
nflBiBB|BBB5SBBBBB||flRBiiJ«UU#<^tl'#«;jjBBBBBBB 



— e ■ , . * i i, i 

.-oo.*-/i#yEqS|ir- 

— ■.'>A^UBR•J!^^e- 

,* A#W»f!fiBSa' 



i88BSBSSBBBSBBB8BBShBBii»«ffUUU^MWWU 

■ ■■■RRlieBBiBBIBBBBBBB 

BBIBBBIBBBBBBBSBBB 



IBBB 



ie§BBBBBSBBSiBSS8BSBBB8BBB{Sa[yUy^ 



1.273YBBBBBBB 
IBBBBB 



BRBBR 



1.000 



1.610 



2.238 



2.966 



3.6<»<. 



A. 322 



5.000 



«:rAi P FArinR*; ■ X-AXlSi E+00 Y-AXISt E + 00 Z-AXIS« E + 00 . ^. 

15^2?,/*^^° 0( -9)/ 1.003 (-9), 2.000(-9), 3.000 (-9 ), ^-ROO -9 

Z5-Z9» 5.000i-9)i 6.000 (-9), 7.000(-9), 8.000(-9 ), 9.000(-9) 

'ROGRAM -TEST- REiDY FOR PUT 



FINIS ■ .T.* $ 

ROGRAM -TEST- TERMINATED 
3R 



STOP PROGRAM 



6 



3^ 
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APPENDIX: ENRICHMENT 
One dimensional enrichment is applicable only to data sets whose 
points describe an ordered sequence of line segments, as follows: 



^ ^ ^ : Line 1 



ERIC 



'^i+l' ^i+1* % 



Line 2 



'^j+l' ^j+l* ^c 



c 



One dimensional enrichment assumes that neighboring pairs of points whose 
z values are identical are part of the same line, hence it generates a 
sequence of equi-spaced points to fill in the gap between them. Note that 
the z values within a line must be absolutely identical; not merely similar, 
hence they must be set by means of assignment statements and not as the 
result of floating point computations. 

The enrichment algorithm operates as follows. Let 

^i " (^i^yi'^i>' ^+1 " ^^i+r^i+r^i+i^' ^+2" K+2'yi+2'^i+2^ 

97 
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be three sequential points in the binary data file where = z^^^^ = z^^j* 

The enrichment algorithm adds the points 

. Qj " + jAx^, + jAy^, z) for j = 1, . . . , n^ 

between and P^^^.* points 

% " ^^i+1 ^^2' ^i+l ^^2* z) foi^ k « 1, ... , n^ 
between P^_^^ and P^_^2 where 

^ " = ^i+l = ^+2' 

= Max[ I (x^^^-x^) | /XRICH, | iy^+j-y^) \ /YRICH] , 

n2 = Max[|(x^_^2-Vl^l^^^^"' I (y^+2"yi+i^ 1^^^"^ » 















Ax2 = 






Ay2 = 


^^1+2 





and XRICH, YRICH are values specified by the user to control the spacing 
of enrichment points (the number of points is inversely proportional to 
the magnitudes of XRICH and YRICH). Both XRICH and YRICH must be specified 
because the x and y axes of the plot may have radically different scales. 
If z^ « ^i+l ^i+1 ^ ^i+2' ^^^^ enrichment assumes that P^^^. t^^i'^^t^s 
a line and ^^^^2 ^^^^^^ ^ line, hence enrichment takes place only 
between P^ and '^^^y If z^ ^ ^±-^-1 ^± " ^i+2* ^^^^ enrichm^^nt takes 
place because no neighboring z values are identical. 



38 
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For example, if XR1CH«10 and YRICH=1 and the file of x, y, z coordinates 



is as follows: 

Sequence No. x Coord y Coord z Coord 

1 8 5 A2.0 

2 30 7 A2.0 

3 35 A 42. 0 

A 50 8 42.001 

5 20 3 42.0 



then PRNT3D adds a series of linearly interpolated points to the file 
between points 1 and 2 and 2 and 3, but not between 3 and 4 or 4 and 5. 
The number of points, and the Ax and Ay Increments for the points added to 
lines 1-2 and 2-3, are as follows: 

nj = max ( 1 30-81/10, |7-5|/1) -= 2 
n^ = max ( 1 35-30 1 /lO, |4-7|/1) = 3 
Ax J - (30-8) /2 = 11 
Axj = (35-30)/3 = 1.66 
Ay J = (7-5) /2 - 1 
Ay^ = (4-7) /3 - -1 

The set of points added to line 1-2 are: 

(Xj+Axj, yj+Ayj, z) " (8+11, 5+1, 42) - (19, 6, 42) 
(Xj+2Ax,, yi+2Ay,, z) - (8+22, 5+2, 42) = (30, 7, 42) 

and the points added to line 2-3 are: ^ 
(Xj+Axj, yj+Ay^, z) - (3C>fl.66, 7-1, 42) - (31.66. 6, 42)^ 
(X2+2AX2, y2+2Ay2, z) - (3CH-3.32, 7-2, 42) - (33.32, 5. 42) 
(Xj+SAXj, y2+3Ay2. 2) - (30+4.48, 7-3, 42) - (34.48, 4, 42) 



Two-dimensional enrichment is applicable only to data files which 
are organized in a special way; namely, the data points in a file must 
form a set of triangles in the x,y plane. Consider the following region 
.which is divided into four contiguous triangles labeled T^,...,T^ whose 
vertex points are labeled P. 




Each triangle is defined by three sets of x,y,z coordinates and the data 
file must be organized as follows: 

P^: x^, y^. 

^2' ^^2' ^2' ^2 • ^1 

P5: X5, yy Z5 



^2' ^2' ^2'. ^2 

P3: X3, yy Z3 : T2 

P5: X5, y^, 



p^* '^3» y^' ^3 

^4- ''A' ^A' • ^3 

P5: x^, y^, 



30 



Note that each point occurs at least twice in the data file. Enrichment 
fits a triangular plane through the z coordinates, of all three points in 
each triangle. The distances between points in this plane are controlled 
by the same parameters XRICH and YRICH as are used in one-dimensional 
enrichment. This form of enrichment or interpolation is particularly 
suitable for displaying the solutions to finite element formulations. 
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